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1.  Introduction 


1 ,3,5-trinitropcrhydro-l  ,3,5-triazinc  (RDX)  and  octahydro-1, 3,5,' 7-tetranitro-l, 3,5,7- 
tetrazocine  (HMX)  are  standard  military  explosives.  RDX  and  HMX  have  deto¬ 
nation  velocities  of  more  than  8.8  km/s,  Chapman- Jouguet  (CJ)  pressures  above 
33  GPa,  and  total  energies  of  detonation  above  10  kJ/mL.  The  performance  of  a 
given  propellant  is  largely  dominated  by  its  CJ  pressure.  It  is  therefore  desirable  to 
find  explosives  with  comparable  detonation  characteristics,  i.e.,  shock  velocity  and 
total  energy  of  detonation,  to  RDX  or  HMX  while  exhibiting  increased  ability  to 
propel  loads.  This  enables  less  propellant  or  higher  final  velocities  in  applications, 
but  the  usefulness  is  curtailed  by  the  synthetic  accessibility,  which  may  result  in 
higher  costs. 

Inverse  molecular  design  (IMD)  is  a  method  for  materials  optimization  well  suited 
for  the  exploration  of  energetic  materials.  Of  late,  IMD  beyond  the  purview  of  drug 
design  has  enjoyed  increasing  popularity.  Examples  can  be  found  in  protein  de¬ 
sign1^  or  high-hyperpolarizability  materials.5'7  The  design  problem  is  complicated 
by  the  vastness  of  possible  chemicals,  termed  chemical  space.  This  space  can  be 
viewed  as  combinatorially  complex,  e.g.,  20s  ~  2.6  •  1010  octa-peptides  of  the  nat¬ 
urally  occurring  amino  acids  exist  alone.  As  a  consequence,  a  variety  of  methods 
have  been  developed  for  the  discrete  optimization  in  chemical  subspaces.8'15  The 
continuous  optimization  of  chemicals  used  in  the  linear  combination  of  atomic  po¬ 
tentials  (LCAP)  method1617  and  the  variation-of-particles  density-functional-theo¬ 
retical  (VP-DFT)  method  introduces  an  important  concept  for  dealing  with  the  in¬ 
herent  roughness  of  chemical  compound  space.18-20  LCAP  and  VP-DFT  interpolate 
continuously  between  the  Hamiltonians  of  various  chemical  species.  Furthermore, 
recently  an  investigation  into  the  reasons  why  chemical  optimization  is,  relatively 
speaking,  “easy”  utilized  probability  distributions  and  expectations  on  the  control 
variables  to  arrive  at  its  conclusion.21-22  Recently,  it  was  discovered  that  indeed  all 
properties  can  be  optimized  efficiently  only  using  the  charge  density  with  as-yet  un¬ 
known  functionals.23  All  previous  optimizations  were  executed  using  at  most  single 
constraints.  In  this  contribution,  we  introduce  a  flexible  method  for  multiple  non¬ 
linear  (inequality)  constraints  applied  to  a  combinatorially  complex  search  space  of 
energetic  materials. 
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2.  Computational  Chemistry  Protocol 


All  quantum-mechanical  computations  were  performed  using  Gaussian  09.24  All 
geometries  were  preoptimized  with  B3LYP/3-21G  under  constraints  as  described 
elsewhere25  followed  by  a  geometry  optimization  using  B3LYP/3-21G*  (see  Listing 
1  in  the  Appendix  for  specifics). 

Before  further  analysis,  the  conformational  space  was  explored  to  ensure  low-energy 
species.25  Natural-bond-order  analysis  was  used  to  determine  connectivity  and  cre¬ 
ate  SMILES  (simplified  molecular-input  line-entry  system  identifiers)  for  input  into 
SYLVIA  for  synthetic  viability  analysis.26  The  heats  of  formation  and  density  were 
computed  following  Byrd.27  These  served  as  input  to  Cheetah28  for  energetic  char¬ 
acterization.  Listing  2  summarizes  these  steps. 

3.  Optimization  Methodology 

Our  optimization  problem  is  formulated  as  a  constrained  maximization: 


max  P(x) 
xeccs 


s.t.  :  TED(x)  -  9.75  >  0 
SV(x)  -  9  >  0 
5  —  SA(x )  >  0, 


(1) 


where  TED(x)  is  the  total  energy  of  detonation  (TED)  of  compound  x  from  the 
chosen  chemical  subspace  (CCS)  of  chemical  compound  space  discussed  later, 
S V (x)  is  the  shock  velocity,  P(x)  is  the  CJ  pressure,  and  SA(x )  is  the  synthetic  ac¬ 
cessibility.  We  reformulate  problem  1  via  nonnegative  Lagrange  multipliers  A  £ 
for  the  3  constraints  to  the  augmented  Lagrangian  function  £(x,A)  :=  P(x)  — 
A C(x)  as  a  constrained  min-max  problem, 


max  min  P(x )  —  A TC(x), 

xeCCS  AeR3 


(2) 


where  C(x)  is  the  vector  of  constraint  violations,  i.e.,  ??(9.75  —  TED(x)),  r](9  — 
SV (x)) ,  7j(S A(x)  —  5),  and  is  the  violation  function: 


y,  y  >  0 

0,  otherwise 
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We  solve  problem  2  alternatingly  between  the  primal  variable(s)  x  and  the  dual 
viariables  A.  For  fixed  A,  x  is  varied  along  each  substitution  direction  until  no  im¬ 
provement  is  found.  After  all  search  directions  have  been  searched,  the  minimiza¬ 
tion  with  respect  to  A  is  performed  using  the  dual  function  P(d) (A)  =  maxIgCCS  jC(x,  A) 
and  its  approximation 

P(d)W=  max ■  C(x,\),  (3) 

where  CCS'  is  the  subset  of  visited,  i.e.,  already  computed,  molecules.  Due  to 
the  discrete  nature  of  CCS',  P'(d)  is  a  piece-wise  linear  function  of  A  with  general 
derivative  — C(a;(A)),  where  x(X)  :=  arg  maxyeccs'  £(y,  A).  Due  to  the  piece- 
wise  linearity,  A  only  changes  meaningfully  when  x(\)  ^  x(Xr).  Thus,  whenever 
C(x )  7^  0,  A  is  updated  cumulatively  by: 

AA  =aC(x*)  max  {A,  0}  (4) 

A  =  inf  {A'  e  R+\£(x,  X*  +  A 'C{x*))  >  £(x*,  X*  +  \’C(x*))}  (5) 

where  a  >  1,  x*  is  the  currently  active  molecule,  and  A*  are  the  currently  active 

Lagrange  multipliers. 

Choosing  the  update  of  A  along  the  constraint  violation  direction  C(x*)  results  in  a 
steepest  descent  update  with  a  conservative  step  size.  The  step  size  is  chosen  such 
that  either  x*  violates  constraints  the  least  as  well  as  maximizing  P  or  that  there 
exists  x'  e  CCS'  with  an  improved  £{x' ,  A  +  AA)  >  £(x*,  X  +  AA). 

The  chemical  subspace  that  we  investigate  is  depicted  in  Fig.  1.  Each  X  on  1  repre¬ 
sents  a  substitution  site  for  which  one  of  10  chemical  groups  is  attached.  The  overall 
search  space  thus  encompasses  1010  molecules.  Each  such  substitution  site  repre¬ 
sents  an  independent  search  direction.14  The  substitutions  are  enumerated  and  thus 
can  be  viewed  as  the  integer  positions  on  a  circle.  The  optimization  proceeds  with 
a  local  line  search  for  each  substitution  site  in  the  prescribed  order  of  sites.  When 
neighboring  substitutions  in  the  current  direction  to  the  current  iterate  are  inferior, 
the  line  search  is  halted  and  the  next  direction  is  searched. 
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Fig.  1  Optimization  framework.  Each  X  may  be  replaced  by  -H,  -C(0)0Me,  -N02,  -C(0)Me, 
-C(N02)2H,  -C(N02)3,  -CH2OH,  -CN,  2,  or  3  for  a  total  of  1010  possible  molecules. 

Since  generally  there  is  no  a  priori  knowledge  of  the  proper  search  order  of  each 
substituent,  the  initial  assignment  is  generally  unsuited  for  smooth  optimization. 
To  mitigate  this  problem,  the  enumeration  of  substituents  is  reassessed  after  each 
full  cycle  of  local  searches  in  each  direction.  For  each  substituent  the  average  arc 
tangent  of  the  Lagrangian  value  is  computed, 


(6) 


which  is  subsequently  used  to  order  the  substituents  around  the  circle  starting  with 
the  lowest  scoring  substituent  and  placing  substituents  altematingly  to  the  left  and 
right  in  ascending  order.  The  resulting  order  produces  a  smooth  ordering  that  is 
monotonically  increasing  until  the  maximum  is  reached,  and  then  monotonically 
decreasing  until  the  minimum  is  reached,  or  vice  versa  depending  on  which  part  of 
the  circle  one  starts  and  in  which  direction  one  goes.  The  arc  tangent  is  used  because 
some  computations  may  fail  due  to  convergence  issues  in  geometries  or  general 
instability  of  the  molecule.  In  such  a  case,  a  molecule  is  mapped  to  a  Lagrangian 
value  of  —  oo.  The  general  program  flow  is  visualized  in  Fig.  2. 
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Fig.  2  Flowchart  of  algorithm 


4.  Results  and  Discussion 

If  TED  is  optimized  and  only  SA  is  constrained,  then  nitration  of  all  possible  po¬ 
sitions  occurs.  Steric  hindrances  select  these  positions  when  the  maximum  number 
of  nitro-groups  has  been  reached. 

A  different  picture  emerges  under  full  constraints  optimizing  for  CJ  pressure.  In 
this  case,  the  optimization  still  primarily  nitrates,  but  CJ  pressure  is  attenuated  by 
minimizing  steric  interference  leading  to  exceptionally  high  CJ-pressure  for  mod¬ 
erate  synthetic  effort.  The  final  compound  4  (see  Fig.  3)  exhibits  a  CJ  pressure 
of  40.791  GPa,  a  shock  velocity  of  9.436  km/s,  a  total  energy  of  detonation  of 
12.447  kJ/mL,  and  a  synthetic  accessibility  score  of  4.570  out  of  10. 
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Fig.  3  Best  candidate  found 


5^  Conclusions 

We  have  developed  an  efficient  algorithm  for  handling  multiple  constraints  in  the 
optimization  of  substitutional  subspaces  of  chemical  compound  space.  We  propose 
novel  energetic  materials  (see  4  in  Fig.  3).  Analysis  of  the  trajectory  suggests  that 
molecules  with  the  general  feature  of  sterically  relaxed  nitro-groups  on  the  ring 
are  particularly  promising  for  experimental  evaluation.  While  increased  TED  and 
shock  velocity  directly  correlate  with  the  number  of  nitro-groups,  the  CJ  pressure 
decreases  with  internal  steric  stress.  All  highly  performing  molecules  include  a  sub¬ 
stitution  on  one  of  the  ring  nitrogens  while  nitrating  the  far  ring  carbons. 
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Listing  1  energy _run.sh  runs  the  constrained  geometry  preoptimization  and  the  follow-up  full 
optimization.  It  returns  the  final  coordinates. 

# / / bin /bash 

#  Execution  script 

#  remove  extra  files 
EXEC=$PWD /  g  0  9  _r  u  n 

filename  = ‘basename  ${1}  .  dat  ‘ 

mkdir  — p  $1 

if  [  !  — e  "$l/energy"  —a  !  — e  "$l/failed"  ];  then 

cd  $1 

echo  %chk=opt  .  chk  >  pre.com 
cat  .  ./  header  .  com  »  pre.com 
cat  .  ./$l.zmat  »  pre.com 
echo  %chk=opt  .  chk  >  opt.com 
cat  ../ footer  .  com  »  opt.com 
if  [  !  — e  "pre.log"  ];  then 

NORMALEXEC=  ‘  t  a  i  1  pre.log  I  grep  — o  termination  ‘ 
if  [  !  -z  "  SNORMALEXEC "  -o  !  -e  "opt.  chk"  ];  then 
$EXEC  pre.com  pre.log 
fi 
fi 

$EXEC  opt.com  opt  .  g09_out 

NORMA1  EXET=  ‘tail  opt  .  g09_out  I  grep  — o  Normal  I  awk  ’{  print  \ 
$1  }  ’  ‘ 

if  [  "  $NORMALEXEC "  ==  "Normal"  ];  then 

fgrep  opt  .  g09_out  — e  ’SCF  Done’  I  tail  —1  I  awk  ’{  print  \ 
$5  } ’  >  energy 

awk  — f  .  .  /  logcart  .awk  <  opt  .  g09_out  >  xyz 

python  -/ bin  /  retrieve_zmat_from_xyz  .  py  ../Sl.zmat  xyz  \ 

I  grep  c  >  opt.rconsts 

python  -/ bin  /  retrieve_zmat_from_xyz  .  py  ../Sl.zmat  xyz  \ 

I  grep  d  >  opt . rvars 
.  .  /  retrieve_zmat_g03  opt  .  g09_out  opt 
bzip2  opt  .  g09_out 
bzip2  *.chk 
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else 


cd 

In  — s  $1 /energy  $1 . energy 
In  — s  $l/opt.rvars  $l.rvars 
In  — s  $1  /  opt  .  rconsts  $1  .  rconsts 

echo  —1000000  >  result 
touch  failed 
exit  1 
fi 
fi 


Listing  2  proprty_script.sh  computes  the  properties  and  penalties  of  a  molecule. 

# ! / bin /bash 
#  Execution  script 

EXEC=$PWD /  g  0  9  _r  u  n 

filename  =  ‘ basename  ${1}  .  dat  ‘ 

export  LD_LIBRARY_PATH=$HOME/  slash/lib  :  $LD_LIBRARY_PATH 
export  PATH=$PATH :  $HOME/  pkg  s  /  SYLVIA  /bin 

mkdir  — p  $1 

if  [  !  -e  " $1 / failed "  -a  -e  "$l/xyz"  ];  then 

cd  $1 

if  [  !  — e  "result"  ];  then 

bunzip2  *  .  g09_out  .  bz2 
bunzip2  *  out  .  bz2 
if  [  !  — e  "nbo.out"  ];  then 

cp  .  ./  nbo.com  nbo.com 
tail  — n  +3  xyz  »  nbo.com 
echo  »  nbo.com 
$EXEC  nbo.com  nbo.out 
fi 

if  [  !  — e  "  smile,  smi"  ];  then 

.  .  /  NBOtoSMILE  nbo.out  smile  .  smi 
BLA=‘cat  smile,  smi  ‘ 
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echo  $BLA"  ;  "  >  smile  .  smi 

fi 

if  [  !  — e  " sascore . smi "  ];  then 

sylvia  smile. smi  sascore. smi 

fi 

if  [  !  -e  "penalty"  ];  then 

cat  sascore.  smi  I  awk  — F  ’{print  $3—  5"\n0"}’  \ 

I  sort  — g  I  tail  —1  >  penalty 

fi 

if  [  !  — e  "  heat_of_formation  .  out "  ];  then 

rm  — f  heatstring* 

fi 

../smile2heatstring  >  heatstring.dat 
cat  xyz  I  tail  — n  +3  »  heatstring.dat 
.  ./  master_densheat  .  s  heatstring  .  dat 

MASS=‘obprop  xyz  I  grep  mol_weight  I  awk  ’{print  $2 }  ’  ‘ 
COMPOSITION^  ‘  obprop  xyz  I  grep  ormula  I  awk  ’{print  $2 }  ’  ‘ 
if  [  !  — e  " cheetah . out "  ];  then 

cp  . . / cheeta . top  $1 . chi 

python  .  .  /  build_cheetah_input  .  py  heat_of_f ormation  .  out  \ 
density,  out  smile,  smi  »  $1  .  chi 
cat  ../ cheeta  .  bottom  »  $1  .  chi 
source  ~/pkgs / cheetah7 . 0/ bin / cheetah . sh 
~/pkgs  /  cheetah7 . 0/ bin  /  cheetah  $1  .  chi  >  cheetah. out 
fgrep  cheetah  .  out  — e  ’Total  energy  of  detonation  ’\ 

I  awk  ’{print  — 1*$6  }  ’  >  TED 
fgrep  summary,  out  — e  ’Pressure’  I  awk  ’{print  $3}’\ 

I  tail  —1  >  DetPressure 

fgrep  summary,  out  — e  ’Shock  velocity  ’  I  awk  ’{  print  \ 
$4}’  I  tail  -1  >ShockVel 
cat  TED  I  awk  ’{print  9.75 -$1 " \n0"  }  ’  I  sort  -g\ 

I  tail  —1  »  penalty 

cat  ShockVel  I  awk  ’{print  9.0  —  $1 " \n0"  }  ’  I  sort\ 

— g  I  tail  —1  »  penalty 
cp  DetPressure  result 

fi 
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bzip2  * .  out 

cd 

In  — s  $l/result  $1. result 
In  — s  $l/penalty  $1  .  penalty 
cd  $1 
fi 

bzip2  * .  g09_out 

cd 

else 

exit  —1 
fi 


Listing  3  extract_results.py 

import  sys  ,  io 

def  main  ( filename  ,filenamed): 

with  open  (  filename  )  as  f: 
line  =  f  .  readline  () 
line  =  f  .  readline  () 
r  =  line.rsplitQ 
with  open  (  filenamed  )  as  f: 
line  =  f  .  readline  () 
line  =  f  .  readline  () 
r2  =  line.rsplitQ 

print  float  (r[len(r)  —  1])*  float  (r2  [len(r2)  —  1]) 

if _ name _ -  " _ main _ 

main  (  sys  .  argv  [l],sys.argv[2]) 

Listing  4  build_cheeta_input.py 

import  sys 

sys  .  path  .  append  (  ’  /  usr/people/crinders/slash/lib  ’) 
import  io  ,  openbabel  ,  collections 
NumToName  =  {  6:"C",  7:"N",  1 :  "H"  ,  8 :  "O"  } 
def  is_number  (  s  ) : 

try  : 

f lo  at ( s  ) 
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return  True 
except  ValueError: 
return  False 

def  split_formula  (  formula  ) : 

tokens  =  1  i  s  t  (  formula  ) 

last_was_number  =  is_number  ( tokens  [0] ) 
n  =  0 
w  =  "" 
output  =  "" 

for  i  in  range ( len ( tokens  )) : 

if  last_was_number  : 

if  is_n  umber  (tokens  [i]): 
w  +=  tokens [ i ] 

else  : 

output+="  , "  +w 
w  =  tokens [ i ] 

else  : 

if  is_n  umber  (tokens  [i]): 
output+="  ,  "+w 
w=tokens [ i ] 

else  : 

w+=tokens [ i ] 

last_was_number  =  is_number  ( tokens  [  i  ] ) 

output+=  " , "+w 
return  output 

def  composition_of  (mol ) : 

comp  =  collections,  defaultdict(int) 
for  obatom  in  openbabel  .  OBMolAtomlter  ( mol ) : 
comp  [  obatom  .  GetAtomicNum  ( )  ]  +=  1 

rs  =  "" 

for  k  in  comp,  keys  (): 

rs  +=  NumToName[k]+ "  ,  "  +  str  (comp[k])+  "  , " 
if  len  (  rs  )  >0: 
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rs  =  rs  [0  :  len  (  rs  )  —  1] 

return  rs 


def  main  (  filename  ,  filenamed  ,  compo sition  ) : 
conv  =  openbabel  .  OBConversion  ( ) 
conv  .  SetlnAndOutFormats  ( "  smi "  ,  "  mol2  " ) 
mol  =  openbabel  .OBMol() 

conv  .  Add  Option  ("h"  , openbabel  .OBConversion.  OUTOPTIONS) 
conv  .  ReadFile  (mol ,  composition) 

conv  .  Add  Option  ("h"  , openbabel  .OBConversion.  OUTOPTIONS) 
mol .  AddHydrogens  () 
with  open  (  filename  )  as  f: 
line  =  f  .  readline  () 
line  =  f  .  readline  () 
heat  =  line.rsplitQ 
with  open  (  filenamed  )  as  f: 
line  =  f  .  readline  () 
line  =  f  .  readline  () 
dens  =  line  .  r sp li t  () 

print  "new^form  ,  ...my explosive  , ,  float  (dens  [len  (dens)  —  1]), 

"  ,  float  (heat  [len(heat)—  1])*4. 184,"  ,composition_of  (mol ) 


if _ name _ -  " _ main _ 

main(  sys  .  argv  [1]  ,sys.argv[2]  ,sys.  argv  [3  ] ) 

Listing  5  smile2heatstring 

# ! / bin /bash 


#  NBOtoSMILE  arranges  N3  to  [N-]-[N+]#N 

AZIDES=‘obgrep2  -N  ’  [C] (  - [NX2-] -[NX2  +  ] #[NX1  ])  ’  smile  .  smi  I  head  -1‘ 

#  Nitro  groups  are  always  [N+](=0)[0—] 

CNITROS=  ‘  obgrep2  — N  ’  [C](  —  [NX3]  (  =  [0X1  ] )  [  0X1  —  ])  ’  smile  .  smi  I  head  -1‘ 

ONITROS=  ‘  obgrep2  -N  ’  [0](  —  [NX3]  (  =  [0X1  ] )  [  0X1  —  ])  ’  smile  .  smi  I  head  -1‘ 

NNITROS=  ‘  obgrep2  — N  ’  [N]( -[NX3]  (  =  [0X1  ] )  [0X1-])  ’  smile  .  smi  I  head  -1‘ 

#  NO  bound  to  any  atom 

NITROSO=  ‘  obgrep2  -N  ’  *(  -  [NX2]  (  =  [0X1  ] ) )  ’  smile  .  smi  I  head  -V 
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#  All  SP3  carbons 

CSP3=  ‘  obgrep2  — N  ’[#6^3]’  smile  .  smi  I  head 

#  All  hydrogens 

Hs=‘obgrep2  — N  ’  [H]  ’  smile  .  smi  I  head  — 1 

#  Is  nitrogen  but  doesn  ’  t  have  double  or  triple  bonds 

#  and  is  not  negatively  charged 

NSP3=  ‘  obgrep2  -N  ’  [N&!$  ( [N]  =  *)& !  $  ( [N]  #*  )&!$  (  [N-])]  ’  smile,  smi 
head  — 1‘ 

0SP3=  ‘  obgrep2  — N  ’[0X2]’  smile  .  smi  I  head  — 1‘ 

#  every  carbon  that  is  not  in  CSP3 

CPRIME=  ‘  obgrep2  — N  ’  [#6&!$  ([#6A3])]  ’  smile. smi  I  head  -V 

#  every  nitrogen  that  is  not  in  NSP3 

NPRIME=  ‘  obgrep2  -N  ’  [n  ,  $  (N=*) ,  $  (N#*  )  ,  NX2  ]  ’  smile,  smi  I  head  -1 

#  every  oxygen  that  is  not  in  0SP3 

0PRIME=  ‘  obgrep2  -N  ’  [0&!$  ( [0X2] )  ]  ’  smile  .  smi  I  head  -1‘ 

echo  $CSP3  $Hs  $NSP3  $0SP3  $CPRIME  $NPRIME  $OPRIME  $CNITROS\ 
$NNITROS  $ONITROS  $AZIDES  $NITROSO 

Listing  6  Input  to  DiscreteOpt  defining  the  substitution  framework  and  search  directions. 

ChemGroup  ( 

( 

Z( 


(C, 

-3, 

0 

00  , 

-2, 

000 

0, 

-1, 

000.0) 

(C, 

0, 

1 

44, 

-3, 

000 

0, 

-2, 

000.0) 

(N, 

1, 

1 

44, 

0, 

109 

5, 

-3, 

000.0) 

(C, 

2, 

1 

44, 

1, 

109 

5, 

0, 

60.0) 

(C, 

3, 

1 

44, 

2, 

109 

5, 

1, 

-60.0) 

(N, 

4, 

1 

44, 

3, 

109 

5, 

2, 

60.0) 

) 

ReturnConnector() 
Connector  ( 

( 

(0,1  ,5) 

(0.7,  0.0,  0.0) 
(0.0,  109.5,  0.0) 
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(0.0,  0.0,  120.0) 
(0,0,1) 

(0,0,1) 

(0,0,0) 

0 

) 

( 

(1  ,2,0) 

(0.7,  0.0,  0.0) 
(0.0,  109.5,  0.0) 
(0.0,  0.0,  120.0) 
(0,0,1) 

(0,0,1) 

(0,0,0) 

0 

) 

( 

(3,4,2) 

(0.7,  0.0,  0.0) 
(0.0,  109.5,  0.0) 
(0.0,  0.0,  120.0) 
(0,0,1) 

(0,0,1) 

(0,0,0) 

0 

) 

( 

(4,5  ,3) 

(0.7,  0.0,  0.0) 
(0.0,  109.5,  0.0) 
(0.0,  0.0,  120.0) 
(0,0,1) 

(0,0,1) 

(0,0,0) 

0 

) 
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( 

(0,1  ,5) 

(0.7,  0.0,  0.0) 
(0.0,  109.5,  0.0) 
(0.0,  0.0,  -120.0) 
(0,0,1) 

(0,0,1) 

(0,0,0) 

0 

) 

( 

(1  ,2,0) 

(0.7,  0.0,  0.0) 
(0.0,  109.5,  0.0) 
(0.0,  0.0,  -120.0) 
(0,0,1) 

(0,0,1) 

(0,0,0) 

0 

) 

( 

(3,4,2) 

(0.7,  0.0,  0.0) 
(0.0,  109.5,  0.0) 
(0.0,  0.0,  -120.0) 
(0,0,1) 

(0,0,1) 

(0,0,0) 

0 

) 

( 

(4,5  ,3) 

(0.7,  0.0,  0.0) 
(0.0,  109.5,  0.0) 
(0.0,  0.0,  -120.0) 
(0,0,1) 
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(0,0,1) 

(0,0,0) 

0 

) 

( 

(2,3,1) 

(0.7,  0.0,  0.0) 

(0.0,  109.5,  0.0) 

(0.0,  0.0,  -120.0) 

(0,0,1) 

(0,0,1) 

(0,0,0) 

(240) 

) 

( 

(5  ,0,4) 

(0.7,  0.0,  0.0) 

(0.0,  109.5,  0.0) 

(0.0,  0.0,  -120.0) 

(0,0,1) 

(0,0,1) 

(0,0,0) 

(240) 

) 

) 

allowed_groups  ( 

(1  ,2  ,3  ,4  ,5  ,6  ,7  ,8  ,9)  #  H,  C(=0)0Me,  N02,  Ac,  C(N02)X,  CH20H,  CN,\ 
tetraazaCp  ,  maleic  anhydride 

(1,2, 3, 4, 5, 6, 7, 8, 9) 

(1,2, 3, 4, 5, 6, 7, 8, 9) 

(1,2, 3, 4, 5, 6, 7, 8, 9) 

(1,2, 3, 4, 5, 6, 7, 8, 9) 

(1,2, 3, 4, 5, 6, 7, 8, 9) 

(1,2, 3, 4, 5, 6, 7, 8, 9) 

(1,2, 3, 4, 5, 6, 7, 8, 9) 

(1,2, 3, 4, 5, 6, 7, 8, 9) 
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(1,2, 3, 4, 5, 6, 7, 8, 9) 

) 

) 

(#1 

Z( 

(H,  -3,  0.4,  -2,  0.0,  -1,  0.0) 

) 

ReturnConnector() 

Connector  ( ) 
allowed_groups  () 

) 

(#2 

Z( 


(C, 

-3, 

0 

84, 

-2, 

0.0 

-1, 

000.0) 

(O, 

0, 

1 

3  , 

-3, 

120.0  , 

-2, 

000(180)) 

(O, 

0, 

1 

4  , 

1, 

120.0  , 

-3, 

180.0) 

(C, 

2, 

1 

50  , 

0, 

120.47  , 

1, 

0(180)) 

(H, 

3, 

1 

10, 

2, 

109.47  , 

0, 

180) 

(H, 

3, 

1 

10, 

2, 

109.47  , 

4, 

120) 

(H, 

3, 

1 

10, 

2, 

109.47  , 

4,- 

-120) 

) 

ReturnConnectorQ 
Connector  () 
allowed_groups  ( ) 

) 

(#3 

Z( 


(N, 

-3, 

0.74  , 

-2, 

0.0 

-1, 

000.0) 

(O, 

0, 

1.3  , 

-3, 

120.0  , 

-2, 

000) 

(O, 

0, 

1.3  , 

1, 

120.0  , 

-3, 

180.0) 

) 

ReturnConnector() 
Connector  () 
allowed_groups  ( ) 

) 

( 
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Z(#4 

(C, 

-3, 

0.70  , 

-2, 

000.0  , 

-1, 

000.00) 

(O, 

0, 

1.30  , 

-3, 

120.0  , 

-2, 

000.00) 

(C, 

0, 

1.54  , 

1, 

120.0  , 

-3, 

180.00) 

(H, 

2, 

1.09  , 

0, 

109.5  , 

1, 

180.00) 

(H, 

2, 

1.09  , 

0, 

109.5  , 

3, 

120.00) 

(H, 

2, 

1.09  , 

0, 

109.5  , 

4, 

120.00) 

ReturnConnectorQ 
Connector  () 
allowed_groups  ( ) 

) 

(#5 

Z( 


(C, 

-3, 

0 

74, 

-2, 

000.00, 

-1, 

000) 

(N, 

0, 

1 

10, 

-3, 

109.47  , 

-2, 

180) 

(N, 

0, 

1 

10, 

-3, 

109.47  , 

1, 

120) 

(O, 

1, 

1 

3  , 

0, 

120.0  , 

2, 

160) 

(O, 

1, 

1 

3  , 

0, 

120.0  , 

3, 

180.0) 

(O, 

2, 

1 

3  , 

0, 

120.0  , 

1, 

160) 

(O, 

2, 

1 

3  , 

0, 

120.0  , 

5, 

180.0) 

) 

ReturnConnectorQ 
Connector  ( 

( 

(0  ,  — 3 , 1) 

(0.7,  0.0,  0.0) 

(0.0,  109.5,  0.0) 

(0.0,  0.0,  -120.0) 

(0,0,1) 

(0,0,1) 

(0,0,0) 

0 

) 

) 

allowed_groups  ( ( 1  ,  3  ) )  #  H,  N02 
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Z(#6 


(C, 

-3, 

0.74  , 

-2, 

000.00, 

-1, 

000) 

(H, 

0, 

1.10, 

-3, 

109.47  , 

-2, 

180) 

(H, 

0, 

1.10, 

-3, 

109.47  , 

1, 

120) 

(O, 

0, 

1.40  , 

-3, 

109.47  , 

1,- 

-120) 

(H, 

3, 

1.09  , 

0, 

109.47  , 

1, 

180) 

) 

ReturnConnector() 

Connector  () 
allowed_groups  ( ) 

) 

(#7 

Z( 

(C,  -3,  0.70,  -2,  000.0,  -1,  000) 
(N,  0,  1.36,  -3,  170.0,  -2,  180) 

) 

ReturnConnectorQ 
Connector  () 
allowed_groups  ( ) 


Z(#8 

(C, 

-3, 

0 

70, 

-2, 

000 

00, 

-1, 

000) 

(N, 

0, 

1 

44, 

-3, 

126 

00, 

-2, 

1 7  0(  —  1 80)) 

(N, 

1, 

1 

44, 

0, 

108 

00, 

-3, 

180) 

(N, 

2, 

1 

44, 

1, 

108 

00, 

0, 

000) 

(N, 

3, 

1 

44, 

2, 

108 

00, 

1, 

000) 

(H, 

4, 

1 

44, 

3, 

126 

00, 

2, 

170) 

) 

ReturnConnectorQ 
Connector  () 
allowed_groups  ( ) 

) 

( 
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Z(#9 

(N, 

-3, 

0.7, 

-2, 

000.00, 

-1, 

000) 

(C, 

0, 

1.5  , 

-3, 

126.00, 

-2, 

170) 

(C, 

1, 

1.5  , 

0, 

108.00, 

-3, 

180) 

(C, 

2, 

1.5  , 

1, 

108.00, 

0, 

000) 

(C, 

3, 

1.5  , 

2, 

108.00, 

1, 

000) 

(0, 

1, 

1.4, 

2, 

126.00, 

0, 

180) 

(0, 

4, 

1.4, 

3, 

126.00, 

0, 

180) 

(H, 

2, 

1.1  , 

1, 

126.00, 

0, 

180) 

(H, 

) 

Retv 

Com 

alio 

3, 

1.1  , 

4, 

126.00, 

0, 

180) 

irnConnector  ( ) 
nector () 
wed_groups ( ) 

) 
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List  of  Symbols,  Abbreviations,  and  Acronyms 


CCS 

CJ 

HMX 

IMD 

LCAP 

RDX 

TED 

SMILES 

VP-DFT 


chosen  chemical  subspace 
Chapman-Jouguet 

octahydro- 1 ,3 ,5 ,7 -tetranitro- 1 ,3 ,5 ,7 -tetrazocine 
inverse  molecular  design 
linear  combination  of  atomic  potentials 
1 ,3,5-trinitroperhydro- 1 ,3,5-triazine 
total  energy  of  detonation 

simplified  molecular-input  line-entry  system  identifiers 
variation-of-particles  density-functional-theoretical 


26 


1  DEFENSE  TECHNICAL 
(PDF)  INFORMATION  CTR 

DTIC  OCA 

2  DIRECTOR 

(PDF)  US  ARMY  RESEARCH  LAB 
RDRL  CIO  LL 

IMAL  HRA  MAIL  &  RECORDS  MGMT 

1  GOVT  PRINTG  OFC 
(PDF)  A  MALHOTRA 

1  DIR  USARL 
(PDF)  RDRL  WMM  G 
J  ORLICKI 


27 


Intentionally  left  blank. 


28 


